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figure 1. The test rig situated at Weir Pumps Research Laboratory at Alloa consists of a rigid shaft 
supported by rolling bearings at both ends and a very rigid housing mounted on coil springs, 
(initially, the casing was mounted on air bellows). A pair of simulated balance drum test seals is 
installed in the housing and water is charged from the centre of die casing and discharged from 
both ends. A significant feature of die flow path is the radial inward motion before entry to the 
annular clearance space. Two hydraulic actuators offset spatially at 90° are used to apply forces to 
the middle of the casing. Figure 2 shows a schematic view of the primary test section and figure 3 
shows some details of the rotor and seals. The oscillatory motion of the casing, when excited by 
external forces is modified by the dynamic characteristics of the flow passing through the seals. 

A linearised mathematical model is generally the appropriate starting point to study dynamic forces 
in pump annular seals. Hence the following model was assumed for the seal reaction forces: 



Here K, C, and M are the stiffness, damping, and inertia matrices of the fluid as it passes through 
the clearance space and (x,y) is the instantaneous displacement vector of the seal stator relative to 

the rotor. . . . 

Following the commissioning stage of the test rig, preliminary tests were carried out for moderate 

pressures and rotor speeds and the obtained results were presented in (2). These results showed 
good agreement between theoretical and experimental direct and cross-coupled stiffness 
coefficients. However, experimental direct damping coefficients showed large discrepancies when 
compared with theoretical values. While the theory predicts an increase in damping with increasing 
pressure, the experimental results showed no discernible increase in damping. Cross-coupling 
damping and to a lesser extent fluid inertia could not be estimated accurately. On the other hand, 
measured radial pressure distributions lead to inlet swirl ratios of about 0.8 to 1 .0 as compared to 
the commonly assumed value of 0.5 . 

Now although these preliminary results were seen as very encouraging, the major drawback was 
that the tests could not be extended to higher pressures and higher rotor speeds. The low stiffness 
of the casing suspension system resulted in flow induced vibrations which limited the testing range 
to a rotor speed of about 1600 rpm and a pressure of 4.5 bar. It was clear from the tests that these 
flow induced vibrations were associated with the rotation of the shaft because of their synchronous 
nature and there was strong evidence that they were caused by fluid cross-coupled forces in die 
seals. Using the idea of Massey (3), grooved rings were installed in die axial gap between the back 
shroud of the simulated impeller and the casing with the expectation that the grooves would result 
in a lower fluid swirl at seal entry. They were 25 radial grooves of 45 mm length, 6 mm width and 
1.5 mm depth. Following this simple modification, a significant reduction in die vibration level was 
achieved. In fact, with this configuration the rotor could be run up to die maximum achievable 
speed of 3000 rpm (limited by the DC. motor) at any pressure up to 15.7 bar (limited by the 
delivery pump) with peak displacements not exceeding 25 pm. It is well known that a reduction of 
swirl at seal entry results in a reduction of the fluid cross-coupled forces. It is thus obvious that the 
vibrations experienced without the grooves are a direct result of fluid cross-coupling in the seals. 
Since the test rig design enabled different test rings to be tested, it was possible to investigate the 
effect of various groove designs on the dynamic behaviour of the system. 
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values, differentiating the displacement values, and taking the average. The processed data was 
then used to produce estimates of all die coefficients in the linear model and the appropriate 
standard errors. A data acquisition and reduction flow chart is shown in figure 4. More details are 
given by Brown and Ismail (2). 

3. SWIRL MEASUREMENTS 
3.1. Introduction 

The equilibrium equation in the radial direction, in die axial gap h made the simulated impeller and 
the casing (figure 3) is given by: 

y jK — _ j-jfe ^nw fl 

r dr r p dr p h ph 

Neglecting the radial shear stresses Tp^i and and neglecting the change in radial velocity, 
the above equation becomes: 

prto J 

dr r 

where co is the fluid angular velocity at radius r. Stepanoff (4) gives results of the integration of the 
above equation for a velocity distribution given by an equation of the form; 


to = Cr" 

where C and m are constants, m » 2 for a free vortex and m = 0 for a forced vortex. The well 
known free and forced vortex are special cases of vortex motion. 

The fluid angular velocity at seal inlet is obtained from the measured radial pressure distribution. 
For instance, in the case of a forced vortex we have: 

p{r) = j dp = pco 2 J rdr 

If the pressure is known at a radius rj and at seal inlet (r=R), integration of the above equations 
yields the fluid angular velocity at seal inlet 


co® 


2 /?(fl)-p(Q 

p 


3.2. Radial pressure distribution and fluid swirl at inlet : Results and discussion 


The study of the radial pressure distribution of the flow in the gap between the back of die 
simulated impeller and the casing involved the measurement of die differential pressure at five 
radial positions (figure 3) for various flow rates and shaft speeds. A non dimensional flow 
coefficient Cq defined as : 


was used to analyse the results. 
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If the shaft speed is increased while the flow rate is kept constant, die change in stiffhess in the 
above expression is mainly influenced by the second term which contains P 2 > ^ ,s ^ 
inertia coefficient. However, published experimental work has shown that this term is generally 
over-estimated by the theory. The plotted data also shows that the > theory under-estimates the 
stiffness by about 20 to 25%. A direct measurement (static test ) of stiffness for no rotation is 
shown in Figure 6.c . It is clear that the values are comparable with the identified stiffness. 

4.2*2. Cross-coupled stiffness coefficients 

The theoretical and experimental cross-coupled stiffness coefficients are shown in figure 7. As 
predicted by the theory >7 experimental coefficients increase linearly with the shaft speed. The 
sign of the coefficients is also detected accurately. There is good agreement behveen fceory and 
experiment at low axial Reynolds number, figure 7.a, the theory overpredicting the coefficients by 
about 10%. But when the axial Reynolds number increases the agreement deteriorates. The theory 
predicts larger coefficients than the experimental ones especially at high shaft speeds, figure 7.b. 
For this case, the predicted coefficients are about twice as large as the experimental ones. In order 
to explain this discrepancy some of the theoretical assumptions were reviewed. 

The theoretical stiffhess expression is based on a model that assumes tiiat (i) the seal is c ™C€»tnc, 
<ii) the rotating and stationary walls of the seal have the same surface roughness and (in) the 
circumferential fluid velocity starts from an arbitrary initial swirl and asymptotically approaches 
0.5RO as it proceeds axially along the seal. In (6), Childs and Kim used a combmed analytical- 
computational method to predict the fluid circumferential velocity in seals which may have 
different surface roughness treatments on the stator or rotor seal elements. They predicted that e 
circumferential velocity solution in a seal with the same rotor and stator roughness converges 
asymptotically towards 0.5Rft irrespective of whether both wall surfaces are smooth or rough. 
However, they predicted an asymptote less than 0.5Rft if the stator is rougher than die rotor and 
greater than 0.5RQ if the rotor is rougher than the stator. Further evidence that the fluid 
circumferential velocity inside a seal can be different from 0.5R ft was provided by Addlesee et al 
(5). Addlesee et al used a specially designed test rig to study the distribution of fluid velocity at 
different axial positions in a rotating annulus utilising an LDA technique. The test results reported 
in (5) suggest that the fluid circumferential velocity is strongly dependent on the inlet swirl. 
Specifically, strong inlet swirl gave high circumferential velocities (greater than 0.5 Rft) and weak 
inlet swirl gave low circumferential velocities (less than 0.5 RQ). On the other hand, the fluid 
cross-coupled stiffhess varies proportionally with the fluid circumferential velocity. Since die mtet 
swirl in this work can be considered as weak and that the stator surface roughness is greater than 
die rotor surface roughness (1.5 pm and 0.5 pm), it can be argued that the discrepancy between 
theoretical and experimental cross^oupled stiffness coefficients is the result of a combmed effect 
of a rougher stator and a reduced fluid circumferential velocity. This view was supported by a 
good correlation between predicted and measured cross-coupled stiffhess coefficients when smaller 
values of fluid circumferential velocities were assumed. Larger discrepancies between predicted 
and measured cross-coupled stiffhess coefficients at the higher speeds are consistent with flow 
measurements which showed that higher speeds lead to smaller inlet swirl ratios, resulting m 
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smaller values for the fluid circumferential velocity. An other factor that contributes to the 
reduction in fluid circumferential velocity is the axial velocity of the flow. When the fluid enters the 
seal with an initial swirl, the axial pressure drop imposes an axial velocity component on the flow 
which then proceeds in the seal with a combined peripheral and axial motion. But as the axial 
velocity increases, the passage time T becomes smaller and die fluid motion becomes dominated by 
the axial velocity. As a result, a fluid particle entering the seal with some initial swirl may exit it 
with little peripheral motion or more specifically with a weaker peripheral motion than in the case 
of a larger passage time. 

42 J Direct damping coefficients 

The experimental and theoretical direct damping coefficients are shown in figure 8 for R a values of 
6950 and 12 120. As predicted by the theory, the direct damping does not vary with shaft speed. 
The theoretical and experimental coefficients coincide at low axial Reynolds number but disagree 
at high axial Reynolds number. The experimental data indicates that the damping is almost 
insensitive to changes in Reynolds number, which is in contradiction with the theory. It is 
worthwhile noting that similar discrepancies between theoretical and experimental direct damping 
coefficients were also reported by Childs et al in (6), (7) and (8). In (6), Childs and Kim reported 
that measured damping coefficients were smaller than (he theoretical ones at low axial Reynolds 
number but as the Reynolds number increased the agreement between theory and experiment 
improved. 

42.4 Overall dynamic behaviour 

The linear modelling of the complete system was tested by comparing the simulated exciting forces 
from the linear model with the experimentally measured forces flora the two hydraulic jacks. As 
can be seen from figure 9 the comparison demonstrates that an overall RMS error of 16% gives a 
very good reproduction of the forces acting on the casing. Thus the linear model is shown to be 
capable of predicting vibration response accurately 

5. CONCLUSIONS 

Extensive testing with the grooved plates allowed a constructive critical review of the theoretical 
model to be made. Reassessment of the theoretical assumptions was made possible by carrying out 
tests covering the span of die flow rate range permitted by the supply pump. Testing at different 
flow rate values allowed some of the flow underlying phenomena to be unravelled. The obtained 
results are summarised in the following section: 

• The inlet swirl ratio co /ft increases with the flow rate and decreases with the shaft speed, the 
sensitivity to flow rate being smaller in the high speed range. 

• Radial pressure measurements in the gap between the back of the simulated impeller and 
the casing show pressure distributions which suggest the occurrence of a forced vortex 
motion, especially at high shaft speeds which is typical of modem centrifugal pumps. 

• The use of radial grooves of 1 .5 mm depth at the entrance to the balance drum resulted in a 
significant reduction in inlet swirl. 

• The measurements show that only the direct stiffness, cross-coupled stiffness and direct 
damping coefficients are important. The cross-coupled damping coefficients were found to 
be insignificant as the least-squares technique gave standard errors which are larger than the 
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coefficients themselves. The inertia terms were identified but the scatter was large. These 
terms were found to lie in the range 35 to 150 Kg as compared to the theoretical value of 87 
Kg. The plotted data showed no identifiable trend. 

• Evidence of the beneficial effect of the grooves was noticed early when the dynamic 
instability experienced with the plain plates was completely suppressed. Later, 
experimental measurements of the dynamic characteristics of the seals showed a significant 
decrease in cross-coupled stiffness with the 1.5 mm deep grooved plates. 

• The cross-coupled stiffness coefficients were fairly well predicted by the theory at low flow 
rates and small groove depth but the agreement between theory and experiment deteriorated 
with the increase of flow rate and groove depth. The theoretical model was reassessed by 
assuming that the fluid circumferential velocity in the seal is much smaller than half shaft 
speed. 

• The direct stiffness coefficients were found to be under-predicted by the theory by about 20 
to 30% in general and were found to be insensitive to swirl reductions. The measured 
coefficients are supported by those obtained from direct static tests. 

• The measured direct damping coefficients were found to be insensitive to shaft speed, flow 
rate and swirl reductions. This is in serious discrepancy with the theory which predicts an 
increase in damping with increasing flow rate. However, impact tests carried out at the 
highest flow rate showed extremely good agreement with the identified damping. 

Finally, it is to be noted that all identified coefficients except the cross-coupled damping terms 
were identified with small standard errors. The calculated root-mean square errors (RMSE) 
between the measured forces and those obtained by fitting the experimental data to a ten-coefficient 
model were of the order of 8 to 20%. The RMSE is defined as 


RMSE(f) = 100 x 


where / is the measured force vector, / the force vector predicted by the model, c f the standard 
deviation of the measured force and N is the number of sample points in the data records. The 
RMSE provides a measure of the goodness of fit between the experimental data and the 
theoretical model. 
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NOTATION 

seal length, seal radial clearance, axial gap between wall and rotating disk 
radial position, seal radius, simulated impeller radius 
components of the external force vector in x and y directions 
seal restoring forces in x and y directions 
average fluid axial velocity, average fluid radial velocity 

fluid circumferential velocity, fluid circumferential velocity at seal entrance 
time (L/V) 

stiffness, damping and inertia coefficients 
pressure, flow rate, flow coefficient 
average axial Reynolds number, 2pVc/p 
average circumferential Reynolds number, pftRc/p 
friction factor, o » XL/c 
fluid Viscosity, fluid density 
pressure drop across seal, inlet loss coefficient 
shaft angular speed, fluid angular velocity, 

T rwall * ^rdisk radial wall shear stress, radial shear stress on rotating disk 


L, c,h 

r , R , R d , 

Fx.Fy 

R x ,R y 

V, V 

V C’ V 0 

T 

K, C, m 
P.Q.C q 

Ra 
Rr 
X , o 

M. P 

ft ,G> 
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DISCUSSION 


Professor Childs raised a comment concerning the acceleration terms in the equations of 
motions. The equations of motion for the stator mass M s can be written as 


Fy-M.y. 


Kja Kxj fAx) C** C X y fAxl ® J Air 1 

KfX | Ay J ^ Cy X Cyy ^ ^ Ay J 0 My y | Av j 


Here, (/; ,F f ) are the measured components of the input forces, (x t ,y t ) are the measured 
components of the stator acceleration and (Ax , Ay) are defined bellow: 

Ax = x, - x r 
Ay =y,-y. 


where the subscripts s and r identify the stator and rotor respectively. We measure the 
absolute stator displacements (x, ,y t ) and the relative displacements (Ax, Ay). The measured 
displacements show that when the synchronous noise is removed from the data we have 

Ax = x. 

Ay * y, 


and hence the use of absolute accelerations is justified. The above condition is implied by 
the rigid rotor design used. Should this condition not be satisfied one can use a numerical 
procedure to obtain the relative accelerations from the measured relative displacements 
(Ax, Ay). In fact this procedure is currently being adopted for the high pressure tests. 
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Figure I Test rig layout; instrumentation and hydraulic shakers not shown 
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Bearing homing 


Discharge chamber 



Detalb of shakers attachment 


Figure 2 Schematic view of primary test section 
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1-10 Static pressure tappings 

Cl Casing 
C2 Simulated impeller 
C3 Rotor 

C4 Plate for axial gap control 
C5 Seal carrier 
C6 Bronze sleeve 
C7 anti-swirl plate 


Figure 3 Cross-section of test seals 
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Figure 4 Data acquisition and reduction flow chart 
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Figure 5 Inlet swirl ratio versus flow coefficient C q 



Rotor speed (rpm) 


Figure 6.a Direct stiffness coefficients at R a =6950 



Rotor speed (rpm) 


Figure 6.b Direct stiffness coefficients at R a =12120 
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Koverage=1 1 .2 MN/ m 
• x 4 y 



Figure 6.c Static stiffness of system at R a =12100 and N=0 rpm 



Rotor speed (rpm) 


Figure 7.a Cross-coupled stiffness coefficients at R a =6950 



Rotor speed (rpm) 

Figure 7.b Cross-coupled stiffness coefficients at R a =12120 
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Ra=6950 Groove depth-1 .5 mm 



Rotor speed (rpm) 


Figure 8.a Direct damping coefficients at R a =6950 


Ra = 12120 Groove depth=1 .5 mm 



Rotor speed (rpm) 


Figure 8.b Direct damping coefficients at R a = 12120 



time (s) 


Figure 9 Comparison of measured and modelled forces 
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